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We show that recent available lattice data for the renormalized Polyakov loop above the decon- 
finement phase transition exhibit unequivocal inverse power temperature corrections driven by a 
dimension 2 gluon condensate. This simple ansatz provides a good overall description of the data 
throughout the deconfinement phase until near the critical temperature with just two parameters. 
One of the parameters is consistent with perturbation theory while a second, non perturbative, 
parameter provides a numerical value of the condensate which is close to existing zero and finite 
temperature determinations. 
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The Polyakov loop plays a relevant theoretical role in QCD at finite temperature. It represents the propagator of 
a static test quark and therefore it is crucial in the understanding of the confinement-deconfinement crossover. In 
P, Q it was related to a heavy quark free energy so its vanishing in quenched QCD signals the confinement phase. 
As noted by 't Hooft |3j , gluodynamics at finite temperature formulated using the imaginary time formalism, has an 
7-H ■ extra discrete global symmetry, in addition to usual gauge invariance. This symmetry is spontaneously broken above 
' the deconfinement phase transition 5|, The Polyakov loop, L(T), is a natural order parameter for such phase 
transition; under periodic gauge transformations L is an invariant object but under a 't Hooft transformation it picks 
up a factor which is an element of the center of the gauge group. Effective field theories for the Polyakov loop have 
been proposed in [(|. (For a comprehensive review see e.g. Ref. 0). 

The smooth Wilson loops, and in particular the Polyakov loop, are composite operators. Their perturbative renor- 
' malizability was discussed in pi. IsL [TfJ. ITT| . finding the remarkable result that they are multiplicatively renormalizable, 
without mixing with other operators. Soon afterwards, the perturbative evaluation of the Polyakov loop was ad- 
O ^ dressed by Gava and Jengo [12( within dimensional regularization, to next-to-leading order (NLO). After including 
finite temperature vacuum polarization effects through Debye mass insertion, the leading order term turns out to be 
i-^H C(<? 3 ) instead of the naively expected 0(g 2 ). Their result implies that at high enough temperatures the renormalized 
Polyakov loop should approach unity from above, a consequence of the non trivial factor introduced by the renormal- 
ization. (The expectation value of the bare Polyakov loop vanishes in the continuum limit in any phase.) Not much 
progress has been achieved after this early result. At pres ent there are no perturbative calculations of the expectation 
value of the Polyakov loop beyond NLO. As noted in [121 a direct calculation would have to confront the proliferation 
of Feynman diagrams due to infrared divergences [13|. A different approach, related to the dimensional reduction 
technique, is discussed below. 

On the non perturbative side, the bare Polyakov loop has often been studied numerically within lattice gauge theory 
calculations, however, a reliable definition and calculation of the renormalized Polyakov loop has been achieved only 
recently. The method introduced in Ref. 0] for quenched QCD obtains the Polyakov loop as a byproduct of the 
heavy quark-antiquark potential at finite temperature, obtained from the correlation between two Polyakov loops at 
different separations. Comparison with the zero temperature potential for small separations allows a quite precise 
determination of the quark selfenergy to be removed and so of the Polyakov loop. The renormalized Polyakov loop 
is larger than unity for temperatures at and above 371, , in agreement with the perturbative expectation. The same 
technique has been applied to two flavor QCD in [Tj|. A direct lattice calculation of the Polyakov loop has also 
been reported in Ref. using a different approach. In this case a single Polyakov loop is used. Comparison of 
data taken at different temperatures allows to determine the renormalization factor to be applied to the bare result. 
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The results of these two approaches agree approximately near the phase transition, but for temperatures above 1.3 T c 
the behaviors turn significantly different. The differences could be due to the effects of finite lattice spacing or to 
ambiguities in the renormalization prescription. 

High temperatures probe kinematical regions which up to the manifest breaking of the Lorentz invariance correspond 
to large Euclidean momenta in the zero temperature quantum field theory. In dimensional regularization in the 
MS scheme one finds that to a given temperature T there corresponds an Euclidean scale /i ~ AttT ^tJi so that 
T c = 270 MeV means /i = 3 GeV. In this regime one expects Operator Product Expansion (OPE) ideas to apply 
and more specifically, at not too high temperatures, condensates and power corrections should play a role. Actually, 
following some older proposals |18|. phenomenological requirements 0, theoretical studies [2(| and lattice analyses 
[2ll I22I I23] ] there has recently been mounting evidence that the lowest condensate order BRST invariant condensate is 
of dimension 2. Such a condensate is generally non-local but in the Landau gauge becomes the local operator (A 2 a ), 
with A^ a the gluon field. Also the (Aq) condensate appears as a parameter in the calculation of the pressure at finite 
temperature |24| . 

In this work we investigate the role of condensates on the expectation value of the Polyakov loop. The Polyakov loop 
is closely related to the thermal expectation value of tr(Ag) (the NLO perturbative result can be obtained in this way) 
and so condensate contributions to this quantity would have immediate impact on the Polyakov loop. Our motivation 
is best exposed by drawing an analogy with the zero temperature quark-antiquark potential in quenched QCD. The 
potential is, of course, closely related to the correlation function of two thermal Wilson lines. The perturbative 
regime of the potential V(r) corresponds to small separations, where the potential is approximately Coulombian. At 
separations of the order of I/Aqcd (there is no other scale in gluodynamics) a linearly confining term develops and 
starts becoming dominant. Both pieces of the potential evolve under the renormalization group at a logarithmically 
slow rate. Therefore, modulo radiative corrections, the dimensionless quantity rV(r) is composed of a flat perturbative 
piece plus a power-like term of the type Aq CD r 2 which is non perturbative. In analogy, at high temperatures, we can 
consider the behavior of the dimensionless quantity (tr(Ao)) /T 2 , also directly related to the correlation function of 
two thermal Wilson lines. The analogous of the scale r in the previous case is the scale 1/T here, and certainly for 
large T the quantity (tr(Ag))/T 2 is perturbative and flat modulo a logarithmic dependence. At lower temperatures 
we contemplate the possibility of non perturbative power- like terms of the type Aq CD /T 2 to develop. As we show 
in this work, such term enters naturally through OPE corrections to the gluon propagator driven by condensates. 
An analysis of available lattice data turns out to display precisely the power-like pattern expected from the previous 
considerations. The pattern is followed in the deconfinement phase from the highest temperatures available down to 
near to the transition where deviations start to show up. 

The paper is organized as follows. In section [H] we discuss perturbative aspects of the Polyakov loop and the 
use of dimensional reduction to attempt the calculation beyond NLO. In section ITTT1 we show that the presence of 
condensates introduce a power-like pattern in the logarithm of the Polyakov loop expectation value. In section HVI we 
analyze the lattice data and show that they are fairly well described as a composition of perturbative plus condensate 
contributions. Finally, in section we summarize our conclusions. 



where ( ) denotes vacuum expectation value, tr is the (fundamental) color trace, and P denotes path ordering. Aq is 
the gluon field in the (Euclidean) time direction, Aq = ^T a ^4o,a, T a being the Hermitian generators of the SU(iV c ) 
Lie algebra in the fundamental representation, with the standard normalization tv(T a Tb) = S a f,/2. 

As a composite operator the Polyakov loop is subject to renormalization. The multiplicative renormalizability of 
the Polyakov loop was established in Refs. ja, EJ fToL lllf in the context of perturbation theory. Gava and Jengo [l2T ] 
addressed the perturbative computation of L(T) in pure gluodynamics. The calculation was carried out to NLO, 
which corresponds to 0(g i ), using dimensional regularization and in the Landau gauge. The result is of course gauge 
invariant. Explicitly, 



II. THE PERTURBATIVE POLYAKOV LOOP 



A. 



Perturbative results 



The (expectation value of the) Polyakov loop is defined as 




(2.1) 



L(T) = I + 




. 9 4 (log^ + J) + 0( 5 5 ). 



(2.2) 
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Here mo is the Debye mass, which controls the screening of chromoelectric modes in the plasma. To one loop |25j 

m D = gT(N c /3 + Nf/Q) 1 ' 2 , (2.3) 

N c being the number of colors and Nf the number of flavors, to account for dynamical quarks. The coupling constant 
g runs with the temperature following the standard renormalization group analysis and one expects (|2.2|) to hold 
for high enough temperature. Remarkably, L(T) turns out to be larger than unity implying that the renormalized 
Polyakov loop is not a unimodular matrix. Note that tud contains a g and so the first non trivial contribution to 
L is 0(g 3 ), due to the infrared structure of the theory, rather than the naively expected 0(g 2 ). Note also that the 
perturbative result (as well as m^) has a well defined large N c limit, with 't Hooft prescription of keeping g 2 N c fixed. 



B. Dimensional reduction 



The result just quoted is rather old yet no higher order computations are presently available. Most efforts in 
perturbative high temperature QCD have been addressed to obtain the pressure and only recently such computations 
have been taken to their highest possible perturbative order j26j, using dimensional reduction ideas |25ll27l l28. 29, 
In order to subsequently include possible contributions from condensates, we will presently reproduce the lowest order 
perturbative result for L(T) using the dimensional reduction approach. In addition this will allow us to discuss 
properties of higher order perturbative contributions to L(T). 

The starting point is the Euclidean QCD action (D M = 8^— igoA^, F^ v = ig^ 1 ^^, D v ], Nf massless fermions) 

£qcd = ^tr(F^) + i> p V + £ gf+g h+ct . ( 2 - 4 ) 

where £ g f+ g h+ c t accounts for gauge fixing and ghost terms as well as the counterterms for renormalization. Next, one 
proceeds to integrate out the fermionic modes and all non stationary gluon modes, which become very heavy at high 
temperature. This results in an effective theory for the remaining stationary (time-independent) gluon modes A, Ax) , 
described by a three dimensional action J d 3 x C 3 (x). To one loop and in the Landau gauge one obtains |l7l l30l l3ll l3a| 

TC 3 (x) = m ^tr(^) + 4^(tr(^)) 2 + ^4(^-^)^(4) 



4tt 2 v v k) " 12tt 2 
tr([A, A } 2 ) + JL^ltr(F§) + T5C 3 (2.5) 



where g([i) is the running coupling constant in the MS scheme (to be used in the Debye mass and in the Polyakov 
loop formula too) 

= 2/3 tog(n/A m ) , 0o = (llJV c /3 - 2AT / /3)/(4 7 r) 2 (2.6) 

and 

' ' 2/3 (log( M /47rT) + lE ) + — ^ (N c + 8N f (log 2 - 1/4)) , 



g%(T) g*fr) ^ 6 ^ > laj 3(4vr) 2 
^T) = m -^o(logW4 7 rT )+7£) + ^ F (-AT c + 8AT.log 2) . (2.7) 

The remainder ££3 contains operators of mass dimension 6 and higher. In addition there are higher loop terms and 
constant (field independent) terms which would be relevant for the pressure. (Note that gE and gM are not to be 
confused with the coupling constants under the same name appearing, e.g., in |33|.') 

At lowest order we will only need the mass term and the kinetic energy term of the chromoelectric field (first and 
fourth terms respectively in Eq. (|2.5Jl ). It will be convenient to work with a rescaled Ag field equal to ff(/i)/5fi(T) 
times the MS A$ field. To all effects, including the Debye mass and the Polyakov loop formula which depends on the 
product gAo, this is equivalent to using the new Aq field together with gE(T) as coupling constant. The latter will 
be denoted g(T) or just g from now on, 

C 3 (x) = ^tr(^) + l t r([A, A ?) +■■■, (2.8) 



1 



9 2 (T) 



2/3 log(T/A E ), 
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with 

For computing the QCD pressure one can use any gauge fixing to integrate the non stationary modes. This is an 
intermediate step to carry out the integration of the remaining modes. Consequently covariant gauges are often used 
as they are computationally simpler. For the Polyakov loop computation the situation is different; static gauges are 
preferred to covariant ones |25| . A static gauge is one in which Aq(x) is brought to be time independent by means of 
a suitable gauge transformation. In such a gauge Eq. (|2.1|1 becomes 

L= ^- (tv e iaA °W T } . (2.10) 

i.e., L depends only on the stationary mode of Aq and so no information is lost on the Polyakov loop operator if the 
non stationary modes are integrated out. Unfortunately, the necessary perturbative computations of e.g. C^{x), are 
only available for covariant gauges. Only in a static gauge the stationary mode A$(x) coincides with the logarithm of 
the Polyakov loop operator. Therefore, in a covariant gauge the effective action of the stationary mode is insufficient 
to recover Polyakov loop expectation values 1 . Nevertheless, as we discuss below, the gauge dependence only affects 
beyond NLO and the two coefficients in Eq. 11'. 21 are reproduced using the formulas in, for instance, [2(|[3i| and the 
method explained in the next subsection. 

Doing a series expansion of L(T) in Eq. (|2.10() one gets 

m = 1 - ^k ma1)) + 4f*w c M)) + • ■ ' • (2 - n) 

tr(^4o) vanish identically while, the other terms of odd order are assumed to vanish due to the QCD conjugation 
symmetry, A^{x) — * —A^(x). The leading contribution is then attached to (ti'(A 2 )}. This quantity has dimensions 
of mass squared and so it would vanish in a perturbative calculation at zero temperature. At finite temperature 
instead it should scale as T 2 modulo slowly varying radiative corrections. Let Doo(k)6 a b denote the momentum space 
propagator for the canonically normalized fields T~ 1 / 2 Ao. a {x), then 

(A 2 o, a ) = {N 2 C - 1)T J —D m {k) . (2.12) 

To lowest order the three dimensional propagator is 

D oo rt (fc) = .2 1 2 , (2-13) 

where the upperscript Pert indicates that it is a perturbative contribution. When this is inserted in (|2.12(l it yields 
(we apply dimensional regularization rules) 

K„r = -(JV C J -1)^ (2.14) 

This result used in Eq. H2.11JI (and using tr(Ap) = Aq a /2) reproduces the perturbative value of L(T) to 0(g 3 ). 

In Fig. ^ we compare the perturbative L(T) in Eq. (|2.2|) with a recent lattice determination of this quantity in 
pure gluodynamics and N c — 3 14]. As we can see, in the high temperature region, T about 6T C , the L(T)-lattice 
is larger than unity, as predicted by the perturbative calculation, moreover the numerical value is also consistent 
with perturbation theory. The agreement quickly deteriorates as the critical temperature is approached from above; 
while the lattice data moves downwards, to eventually displaying a phase transition, the perturbative curve increases 
slightly. As expected, the perturbative result is slowly varying with temperature, the variation coming from logarithmic 
radiative corrections. 



1 Using the stationary mode in Eq. I2.1UI amounts to removing the path ordering operator in the definition of the Polyakov loop, rendering 
it a gauge dependent quantity. 
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FIG. 1: The renormalized Polyakov loop versus the temperature, in gluodynamics. Lattice data from Perturbative LO 

and NLO results are shown for comparison. The curve follows from a fit of the parameter b in Eq. (14.411 . 

C. Higher perturbative orders 

Let us now discuss higher order perturbative contributions to L(T). The renormalizable pieces of the three dimen- 
sional Lagrangian are of the form 



£ren 
3 



■tr([A, M 



m 2 tr(Al) 



Ai(tr(^)) 2 +A 2 tr(^) 



(2.15) 



with Aq ~ T^^Ao, m ~ gT, and Ai ~ A 2 ~ g 4 T. In addition, A = d t - ig 3 A l with A % ~ T~ x l 2 Ai and g 3 ~ T x l 2 g. 
For JV C = 2 or iV c = 3 the A 2 term is redundant and one can set A 2 = 0. The vacuum energy density of this theory, 
f(ga,m, Ai), has been computed to four loops in j2(|, with 173, m and Ai as independent parameters. This allows to 
compute (Aq) and (Aq) by taking derivatives of / with respect to m 2 and Ai respectively, to obtain a perturbative 
estimate of the Polyakov loop. The general structure of the vacuum energy density is as follows (2(J 



/(53,m I A 1 ) = ^^/, fe m 4 -^ 2fc Al 



-k-1 



(2.16) 



where I denotes the number of loops and the coefficients fn~ depend logarithmically on m. Consequently, for the 
quantities in the expansion of L(T) one finds 



y>4 



(tr(A 2 )) 



(tr(4)) 



ff 2 df(g 3 ,m,Xi 
T dm 2 

g 4 df{g 3 ,m, Ai 
T 2 9Ai 



3£ 

E E ^ 

£>1 n=^+2 

^E E 9 n 

l>2 n=l+A 



(2.17) 



As can be seen from these formulas, the first missing contribution to L(T) would be 0(g 7 ) from £ — 5 in the (tr(Ag)) 
term. The lowest contribution from (tr(vlo)) at 5 loops is 0(g 9 ) and that from (tr(Ag)), not available from the 
computation, would first start at 0(<? 9 ) at 3 loops. So in principle, one could extend the perturbative result for L(T) 
to 0(g e ). Unfortunately, the matching relations which connect m, g 3 and Ai to the four dimensional QCD parameters 
are only available in covariant gauges for which the relation l|2.10|l does not apply. In particular, the ratio g(n) / 9b(T) 
used above is gauge dependent at 0(g 2 ) from two loop contributions, this would introduce a gauge dependence at 
0(g 5 ) in L(T). 

On the other hand, the non renormalizable terms ££3 ought to be examined as well to determine to which pertur- 
bative order they start contributing to L(T). The leading such terms are schematically of the type [sTL l3^| 



SC 3 = |jtr([A 



F I 2 

1 -<- fit' J 



) + 



T 3/2 



tr(i? 



9_ 

T 



^A 2 F 2 ). 



(2.18) 
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Using the effective relation Di ~ gT, the first term amounts to an 0(g 4 ) correction to the kinetic energy, so it starts 
contributing at 0(g 7 ) as a correction to the LO. The other terms are effectively of higher order. 

Numerically the terms 0(g 5 ) + 0(g 6 ) computed with the available matching relations do not make a substantial 
contribution as they are qualitatively and also quantitatively similar to those in |12|. Again the radiative nature of 
these perturbative terms produces a rather flat logarithmic dependence with the temperature in sharp contrast with 
the lattice data at not too high temperatures. This reinforces the need of non perturbative effects. 



D. Gaussian ansatz 



It is noteworthy that the contribution from (Aq) starts at 0(g 6 ), and so to 0(g 5 ) Aq obeys a Gaussian distribution. 
That is, to this order one can replace 1)2. 10[) with 



exp 



9 2 (Aj a ) 
AN r T 2 



(2.19) 



and so 



(A), a 



Pert 



JVr 



4tt 



-m D T 



N C (N 2 - 1) 



g 2 T 2 ( log 



~2T 



0(g 3 ) 



(2.20) 



This formula holds also in the unquenched theory, since to this order Nf only appears through the Debye mass. 

The Gaussian ansatz becomes correct 0(g 5 ) at high enough temperature where the theory becomes weakly inter- 
acting due to asymptotic freedom. Also, it becomes exact in the large N c limit as higher order connected expectation 
values are suppressed by powers of 1/N C . Note that A\ a scales as (N 2 — 1) and so L has a well defined limit with 



the standard prescription of keeping 
observed in lattice calculations 341 



g N c finite as N c — » oo. A Gaussian distribution for the Polyakov loop has been 
The Gaussian ansatz is in fact equivalent to expanding the exponential, aver- 
aging over color degrees of freedom and finally invoking the vacuum saturation hypothesis ((^4o fe ) = — ^)^-(A 2 ) k ) 
routinely applied in QCD sum rules at zero temperature. In this line, the Wilson loop was discussed in Ref. [35| by 
using the standard dimension 4 gluon condensate yielding for small contours a term proportional to the area squared 
of the contour. The situation has been revisited in Ref. [2|| in the context of dimension 2 condensates yielding an 
area law for small contours. This agrees with the observation in Ref. that dimension 2 condensates, effectively 
would-be tachyonic gluon masses, provide the short range signature of long range confining forces. 



III. CONDENSATE CONTRIBUTIONS TO THE POLYAKOV LOOP 



As shown in Fig. ^ the perturbative contributions to the Polyakov loop expectation value describe only the region 
of very high temperature. This situation is reminiscent of what happens for the heavy quark-antiquark potential in 
QCD at zero temperature, as a function of the quark-antiquark separation. There, perturbation theory describes 
well the short distance region, where the theory is weakly interacting and standard one-gluon exchange produces 
a Coulomb-like potential. At larger distances confinement sets in and a linear potential must be added to account 
for the lattice data [37j ■ As the potential has dimensions of mass, the Coulomb piece does not need a dimensionful 
coefficient. This makes it allowable in perturbation theory, where Aqcd can only appear through logarithmic radiative 
corrections, as in Eq. (|2.6p . On the other hand, the linear confining piece of the potential requires a dimension two 
coefficient, the string tension, which in pure gluodynamics should be Aq CD times a numerical coefficient. At one loop 
this implies a dependence exp(— 1 / f3og 2 (fi)) , the scale fi being related to the quark-antiquark separation r. While such 
contributions are perfectly possible in QCD, they are clearly beyond any finite order in perturbative QCD and can 
only be attained through suitable resummations of the perturbative series (see e.g. [sHl^l). It * s noteworthy that the 
non perturbative dependence on g is not completely arbitrary, namely, it is such that Aqcd appears raised to positive 
integer powers. This finds a natural explanation from the OPE approach, where the non perturbative contributions 
are driven by condensates of concrete local operators. By dimensional counting, the condensate contributions carry a 
corresponding negative power momentum dependence, so they are subdominant at high momentum as compared to 
the purely perturbative terms but become more important at lower momenta, the lower dimensional operators being 
the dominant ones. In this line the confining piece of the zero temperature heavy quark-antiquark potential has been 
addressed phenomenologically by considering the contribution to the gluon propagator of a dimension two condensate, 
namely, {A 2 ^ in the Landau gauge |l8j |. Just by dimensional counting such term produces a linearly confining term 
in the potential [T^ . 
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In this work we want to investigate the effect of low dimensional condensates on the Polyakov loop expectation 
value. The region of high temperatures is weakly interacting and so ideas inspired on the high momentum region of 
the zero temperature theory might be useful here. As shown above, at high temperatures, the Polyakov loop is closely 
related to the expectation value of Aq in a static gauge. Perturbatively, such quantity necessarily scales as T 2 , but non 
perturbatively a further term proportional to Aq CD is allowed. In order to account for non perturbative contributions 
coming from condensates, we will consider adding to the propagator new phenomenological pieces driven by positive 
mass dimension parameters. Specifically, we consider 



with the non perturbative term 



D 00 (k) = D*t(k) + #o N o° nPcrt (fc) (3-1) 



(k +m z D Y 

Such ansatz parallels those made at zero temperature in the presence of condensates [l8l [l9| . This new piece produces 
a non perturbative contribution to (Aq), namely, 

/ a2 \NonPcrt _ C^c — l)^ m G / o o\ 

\ A 0,a — Z • V" 5- " 3 / 

If we assume that rag is temperature independent up to radiative corrections, the condensate will also be temperature 
independent, modulo these corrections. Equivalently, in terms of the condensate 

o__ ™ I A 2 \NonPcrt 
7-jNonPcrt / r,\ _ 8?F m D /o A\ 

V °° [k) -N2-l T (k 2 +m 2 D ) 2 ' ( ] 

Note that a positive condensate (Aq o ) NonPcrt indicates a would-be tachyonic gluon mass — tUq, as in [l9j . 
Adding the two contributions to (A 2 , a ) in Eq. (|2.19|) . one obtains 

" 2l0gL= 2iV> + 2N C T* (3 - 5) 

The fact that, modulo radiative corrections (including running of the coupling and anomalous dimensions), (Aq a ) Pert 
scales as T 2 while (Aq a ^ NonPcrt i s temperature independent, suggests rewriting the previous formula as ^3 



-21ogL = o + 6(-^J (3.6) 

where the parameters a and b are expected to have only a weak temperature dependence. As advertised the non 
perturbative piece introduces a power-like dependence in the temperature which is not present in the perturbative 
calculation. 



IV. COMPARISON WITH LATTICE DATA 



A. Results in gluodynamics 



A reliable determination of the renormalized Polyakov loop in lattice gauge theory has been undertaken only 
recently in Ref. |14j . for pure gluodynamics and N c = 3. This calculation is, of course, fully non perturbative. These 
authors compute the finite temperature correlation function of a heavy quark-antiquark pair for different separations. 
The two Polyakov loops are multiplicatively renormalized by extracting the (temperature dependent but separation 
independent) quark self energy in such a way that at short distances the standard zero temperature quark-antiquark 
potential is reproduced. At large separations the (squared) renormalized Polyakov loop is then obtained. That is, if 
P x denotes the renormalized Polyakov loop operator located at x, 

(P x P y ) = e- c(T) (pk arc P barc ) = e - F « (r ' T)/T — -> L 2 (T) . (4.1) 
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FIG. 2: The logarithmic dependence of the renormalized Polyakov loop in gluodynamics versus the inverse temperature squared 
in units of the critical temperature. Lattice data from 14] . The fits use Eq. 13.611 with a and b adjustable constants and lattice 
data above 1.03 T c for N T = 4 and N T = 8. Purely perturbative LO and NLO results for Nf = are shown for comparison. 



Motivated by the pattern in Eq. (|3.6[1 . the lattice data for — 21ogL(T) are displayed versus (T c /T) 2 in Fig. El 
As we can see the lattice data follow a nearly straight line. This pattern is clearly distinguishable from the much 
flatter dependence predicted by the perturbative calculation, and unequivocally shows a temperature power correction 
characteristic of a dimension 2 condensate. 

Identification of l|3.6[) with the formula (|3.5|) yields the relations 

1 A 2 -l 2 m D A 2 -l 4 / m D 3\ _ . 
3 2 (^o.a) N ° nPcrt = 2iV c T c 2 6. (4.3) 



A fit of the lattice data of the form 



21ogi = a NLO + 6(^) (4.4) 



T 

with the perturbative value of a to NLO and b as a free constant parameter, yields 



(2.20(6), x2/DOF =/ - 75 < ^ = 4 < (4 5) 

\ 2.14(4), X 1 \l-43, iV r = 8. 1 > 

This corresponds to the following value for the condensate 

a/.a \NonPcrt_ / (0.98 ± 0.02 GeV) 2 , A^ T -4, 
^ ^o,J - |(o.97 ± 0.01 GeV) 2 , iV T = 8 . 

In the fit we include lattice data for temperatures 1.03 T c or above. We use T c /A^ — 1.14(4) an d T c = 

270(2) MeV ^l]- Throughout this section we use the running coupling constant obtained from the beta function 
to three loops and A^ in Eq. I|2.9|l as scale parameter. Assuming that the difference between the two lattice results 
is entirely due to finite cutoff effects, and assuming further that the corresponding leading effect goes as 1/N T , we 
obtain the estimate (0.95(4) GeV) 2 for 3 2 (^o,a) N ° nPcrt in the continuum limit. 

We have also considered a fit of the lattice data with both a and b treated as free constant parameters. This 
produces 

-0.27(5), f 1.81(13), y2/D0F _ (1.07, A T = 4, 

-0.23(1), ° \ 1.72(5), X /UW \0.45, N T = 8 . 

The values of x 2 /DOF are slightly better than the NLO prediction of a. Obviously the identification of a with the 
perturbative result will work better at high temperatures. Using Eq. (|4.2() we obtain for the highest temperature 6 T c 

ft NLO _ _o.22(l) (T = 6T C ), (4.8) 
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in qualitative agreement with the fitted values. Note that for this temperature the non perturbative power correction 
does contribute at the few percent level. For lower temperatures the NLO perturbative result evolves faster than the 
fit suggests. At this level of accuracy one should also take into account logarithmic corrections to the value of the 
condensate and eventually some anomalous dimension correction to the condensate. The present data do not allow a 
clean extraction of such fine details. The average value we get for the condensate with constant a is 

2 /4 2 NNonPcrt_ / (0.89 ± 0.03 GeV) 2 , V T = 4, 
9 {A °> a} ~~ [(0.87 ± 0.02 GeV) 2 , N T = 8 , [ J> 

a little lower than before. The corresponding continuum limit estimate results in g 2 (A^ a ) NonPort — (0.84(6) GeV) 2 . 

We have attempted to determine the coefficient of a possible 1/T 4 correction, appending formula l|3.6[l with a term 
c(T c /T) 4 . When we fit the lattice data for N T = 8, this results in 

6 = 2.18(20), c= -0.04 ±0.24, (4.10) 

with x 2 /DOF = 1.89, where we have considered the perturbative value of a to NLO, and 

a = -0.22(2), 6=1.61(24), c= 0.13 ±0.28, (4.11) 

with x 2 /DOF = 0.42, if we treat a as a free constant. The value of c is compatible with zero in any case, and the 
errors overlap with central values for a and 6 of Eqs. (|4.5|l and l|4.7|l respectively. More accurate data are desirable in 
order to identify contributions from condensates of dimension 4. 

It is noteworthy that a fit to the data completely excludes the existence of a term of the form 1/T in log(L(T)). Such 
term would not have a theoretical basis, as no dimension one condensate exists. However, as noted by the authors of 
|14|. there is a ambiguity in their procedure, which corresponds to adding a constant to the zero temperature quark- 
antiquark potential. Such ambiguity translates into an additive ambiguity in Fq q (r,T) in Eq. (|4.1|l . which would give 
rise a term of the type 1/T in log(T(T)). The absence of such term indicates a preference for the Cornell prescription 
adopted in namely, in Vq q (r) ~ v /r + vi + v 2 r to choose v i = |4^ . 

B. Relation with zero temperature condensates 

Although our determination is based on a static gauge, it is tempting to compare with the zero temperature 
condensate g 2 (A 2 l a ) obtained in the Landau gauge in quenched QCD. There, one obtains from the gluon propagator 
(2.4 ± 0.6 GeV) 2 |2l|, from the symmetric three-gluon vertex (3.6 ± 1.2 GeV) 2 and from the tail of the quark 
propagator (2.1 ± 0.1 GeV) 2 and (3.0 — 3.4 GeV) 2 .231. At zero temperature all Lorentz components are sampled 
suggesting a conversion factor of 4 from g 2 (^4 2 a ) to g 2 (A^ a ), but according to 18], in the Landau gauge the total 
condensate scales as D — 1, D being the Euclidean space dimension, suggesting instead a conversion factor of 3. 
Within the uncertainties of the lattice data as well as the theoretical ambiguities, the agreement is remarkable, as the 
two quenched results refer to different temperatures and gauges. Finite temperature results for the pressure in pure 
gluodynamics yield a value (0.93(7) GeV) 2 for g 2 (^4 2 i. a ) NonPort , indicating an overall coherent picture. 2 

C. Unquenched results 

The renormalized Polyakov loop has also been computed in the unquenched case, using the technique described 
above, in Ref. 0] for two flavor QCD. The lattice data are shown in Fig. [31 and they corresponds to N T = 4. In this 
case, the data fall onto a straight line for temperatures 1.15 T c or above. Closer to the transition temperature the 
data start departing from the pattern l|3.6|) , indicating the need of a richer description as the transition is approached 
from above. A fit to the data above 1.15 T c using a NLO yields 

6 = 2.99(12), 5 2 (A 2 a ) NonPcrt = (0.86 ±0.02 GeV) 2 , (4.12) 

with x 2 /DOF = 1.87. We have used T c /A m = 0.77(9) with T c = 202(4) MeV Q and = 261(31) MeV 0. 
The fit has been done with equal weight to all data points and the value of x 2 quoted corresponds to a representative 
error ±0.05 in 21ogL(T), which similar to that for the quenched case. 



2 This value has been obtained from lattice data shown in figure 2 of Ref. |24j. and also from figure 1 of Ref. l-f'fl . in the temperature 
region used in our fits. 
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FIG. 3: The logarithmic dependence of the renormalized Polyakov loop in unquenched QCD with two flavors versus the inverse 
temperature squared in units of the critical temperature. Lattice data from |15f| . The fits use Eq. 13.61 with a and b adjustable 
constants and data above 1.15 T c . Purely perturbative LO and NLO results for Nf = 2 are shown for comparison. 



A fit with a and b as free parameters gives 

a = -0.31(6), 6 = 2.19(13), <? 2 (^ a ) NonPcrt = (0.73 ± 0.03 GeV) 2 , (4.13) 

with x 2 /DOF = 0.25. As in the quenched case, the value of a is consistent with the perturbative value at high 
temperature 

^nlo _ — 0.35(2) (T = 6T C ). (4.14) 

The lattice data show a departure from the linear pattern for temperatures closer to the transition than 1.15 T c . 
Such departure is not well described by adding new condensates of higher dimension and we have been unable to 
extract a condensate of dimension 4 from the data. We quote here the result of appending a term c(T c /T) 4 in Eq. (|4.4|l . 
The fit of the data above 1.0 T c gives b = 2.44(21) and c = 1.07(19) with x 2 /DOF = 12.8. The coefficients b and c 
are highly correlated. 



D. Further quenched lattice data 



Alternative lattice determinations of the renormalized Polyakov loop in pure gluodynamics have been addressed 
more recently in ^(| . These authors follow a different approach as compared to that in [l4[ • They use single Polyakov 
loops which are multiplicatively renormalized by extraction of the quark selfenergy. The latter is determined by 
isolating the cutoff dependent pieces by comparison of different lattice sizes at the same temperature. Unfortunately 
the results of both approaches differ qualitatively, specially for temperatures above 1.3 T c . This is shown in Fig. 0] 
where the two lattice data sets are compared. 

The origin of the discrepancy between the results obtained with the two approaches is presently not clear, although 
lattice artifacts, in particular finite lattice spacing effects, are not completely excluded in [16j as a possible explanation. 
(Of course, there is also the possibility that after closer scrutiny the two definitions used by the two groups correspond 
really to different renormalized operators.) 

In our view the results in [hH would be the more reliable ones because the method used is technically simpler 
and amenable to tests. Indeed, the authors are able to verify that for small separations of the two Polyakov loops 
the standard zero temperature potential is very accurately reproduced as a function of r for all temperatures. This 
is achieved after a single (temperature dependent) global shift is made, to remove the quark selfenergies; this is 
the quantity c(T) in Eq. 1)4. lfl . The contact between the zero and finite temperature potentials is complete for all 
separations between zero and a T dependent radius r(T) related to the Debye mass, thereby allowing a quite precise 
determination of the counterterm c(T) for each temperature. In addition, as noted above, the calculations are carried 
out for two different lattice sizes, N T = 4 and N T — 8 (and also N T = 16 in 01); an d the results for the renormalized 
Polyakov show very small cutoff dependence, implying that the continuum limit has been reached. 
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FIG. 4: The logarithmic dependence of the renormalized Polyakov loop versus the inverse temperature squared in units of the 
critical temperature. Lattice data from [hH Jld|. The fits use Eq. 13.61 with a and b adjustable constants for the first set of 
data (l4||. and Eq. 14.181 for the second one ;1 (j . 



The method in 16] is technically more difficult to implement (quoting the authors, "In practice, our method is 
not quite so trivial") since it requires comparing different lattice sizes at the same physical temperature. Also the 
subtraction of counterterms is more involved, since, using perturbation theory as guidance, the analogous of c(T) is 
expressed as power series of T with coefficients to be fitted to the bare Polyakov loop data. On the other hand, from 
the point of view of the model proposed in the present work, we expect non perturbative corrections to be negligible 
at the highest temperatures considered in the two lattice calculations and only the data in 01 seem to be consistent 
with perturbation theory 0] at those temperatures. 

The method in 0] renormalizes the logarithm of the bare Polyakov loop by using the scheme 

- log£ barc (T) = / div A r + / rcn + f^N- 1 (4.15) 

where N T is the lattice temporal size, and so N T — A/T, A being the inverse lattice spacing, i.e. the lattice cutoff. 
As said, the data in ^(| deviate from those in ^4[, and in particular, do not follow the pattern l|3.6|l for log(L). Let 
us make a speculation assuming that either the removal of the cutoff dependent pieces has not been complete or that 
after removal of the those pieces, finite renormalization terms of the same type as the subtracted ones remain in the 
renormalized data of [H|. 3 Specifically, let us assume that the data follow the pattern 

-21ogL = a NLO +6(^ + Sa- 1 ^ + Sa + Sa 1 ^. (4.16) 

Actually, we find that the data above 1.3 T c can fairly well be accounted for by using this pattern. This is shown in 
Fig. 2] Remarkably, the central value of the slope b turns out to be close to that found previously with the other set 
of data. However, the best fit has large error bars due to the abundance of parameters available. 

5a = 1.8 ±1.8, b= 1.4 ±2.6, 

<5a_i = -1.0 ±3.8, 6 ai = -0.29 ±0.26, (4.17) 

with x 2 /DOF = 0.0349. 

Similar remarks apply to the fit 

-2logL = a + b(^\ + Sa-x^+5a + 8a x ^r, (4.18) 



3 Of course, one could also ask whether the result in Il4l are not contaminated by finite cutoff effects too, and in particular, whether 
the linear pattern displayed in Fig. |4] is not just the consequence of a huge cutoff effect of the type A 2 /T 2 instead of Aq CD /T 2 as 
proposed in this work. This is unlikely, first because the values of the cutoff A in Il4l are much larger than Aqcd and second, because 
the renormalized results are consistent for different lattice sizes, N T = 4 and N T = 8. 
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although in this case a and 8a cannot be determined independently. This gives 

a + Sa = 1.6 ±1.8, b= 1.3 ±2.6, 

5a_i = -1.4±3.8, 6ax = -0.28 ± 0.26 , (4.19) 

with x 2 /DOF = 0.0350. 

We find encouraging that the value of the condensate approximately agrees using the two different lattice data 
sets. Nevertheless, this speculation is not completely conclusive and an agreement between the results of both lattice 
groups would be needed before further consequences could be extracted. 



V. CONCLUSIONS 



There are two main results of our study. First, when suitably analyzed, the lattice data of the renormalized Polyakov 
loop above the deconfinement phase transition show unequivocally the existence of a non perturbative dimension 2 
condensate. Such contributions have not been considered before but they are in fact dominant and allow to describe 
the data in [l4| down to temperatures as close to the transition as 1.03 T c for pure gluodynamics and 1.15 T c for 
two flavors. Furthermore, the numerical value obtained from the Polyakov loop is quite consistent with the value of 
9 2 (Aq a ^ NonPert extracted from the pressure in gluodynamics. 

We have suggested identifying this condensate with the BRST invariant dimension 2 gluon condensate. Our second 
finding is that, for pure gluodynamics, the numerical value of the condensate (Aq ^N 011 ? 01 ^ defined in a static gauge 
and extracted from Polyakov loop data above the deconfinement transition, is remarkably close to the naive estimate 
(A 2 ^ a }/4, measured at zero temperature and in the Landau gauge. These results pose the theoretical challenge of 
establishing the connection outlined in this paper on a firmer ground. In this light the analogy between the zero 
temperature potential and the Polyakov loop noted in the introduction has been pushed forward in |46| by showing 
that the model in Eq. (|3.1|l predicts a relation between the string tension and the slope of the Polyakov loop that is 
empirically satisfied. 

The simple shape L 2 (T) = e~ a ~ b ( Tc / T ) yields L — > as T — > 0, but does not describe the deconfinement phase 
transition. The closest analogy to such transition would be near the inflexion point of L 2 (T), which takes place at a 
temperature Ti = (26/3) 1 / 2 T c . This T. t would agree with T c for a universal geometrical value b = 3/2, which not far 
from the values obtained in this work from quenched QCD lattice data. Nevertheless, this approximate coincidence 
can only be taken as an estimate since the concrete value of the inflexion point depends on whether L 2 (T) or L(T) 
is used, for instance. It is noteworthy that the same shape can also be obtained within the instanton approach at 
finite temperature along the lines of [471. The relation between instantons and dimension 2 gluon condensates at 
zero temperature was suggested in ,4g and fruitfully exploited in recent lattice simulations to extract, via cooling 
techniques, the infrared behavior of the running coupling constant |49j. In this regard, it might be rather interesting 
to isolate the purely non perturbative instantonic contributions on the lattice and determine whether, after cooling, 
the shape e~( a+b ( T °/ T ) )' 2 extends also below the phase transition. 
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